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Abstract 

We investigate two exclusive processes, e + e~ — > rj c + 7 and e + e~ — > J/ip + J/tp, at the center-of- 
momentum energy y/s = 10.58 GeV within the framework of the nonrelativistic QCD factorization 
43 " approach. A class of relativistic corrections is resummed to all orders in the heavy-quark velocity 

qj v and the corrections are large negative. We further improve the prediction by including available 

£>! 

QCD next-to-leading-order corrections and the interference between the QCD and relativistic cor- 
rections. The prediction for a[e + e~ —> ?? c + t] is about 50 fb. In the case of e + e~ — > J/ip + J/ip the 
standard nonrelativistic QCD prediction for the cross section is negative. As an alternative, the 
vector-meson-dominance approach is employed to compute the photon-fragmentation contribution 
of the process, which gives the cross section ~ 1 fb. This is an indication that the uncalculated 
QCD higher-order corrections may be significant. Our results can be tested against the forthcoming 
data from Belle II and super B factories. 
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I. INTRODUCTION 



The nonrelativistic QCD (NRQCD) factorization approach [1] is a systematic theoretical 
tool to study the production and decay of heavy quarkonia. In this approach, the production 
rate for a heavy quarkonium in a high-energy process is expressed as a linear combination 
of NRQCD long-distance matrix elements (LDMEs) for a heavy-quark-antiquark (QQ n ) 
pair into the physical quarkonium H. The velocity-scaling rules of NRQCD classify the 
numerical importance of each LDME relative to the color-singlet one in powers of v, the 
velocity of the quark and antiquark in the QQ n rest frame. The corresponding short- 
distance factors are calculable perturbatively. While the production of a quarkonium in 
hadron colliders involves various channels with the QQ n pair whose quantum number n 
differs from that of H, an exclusive quarkonium production process in e + e~ annihilation 
such as e + e~ — > J /ip + r\ c or e + e~ — >■ J/ip + J/ip is dominated by the color-singlet process in 
which a produced QQ pair has the same quantum number as that of a final-state quarkonium. 
This greatly reduces theoretical uncertainties. The factorization theorem for some exclusive 
quarkonium production processes in e + e~ annihilation is available [2], although that for the 
inclusive quarkonium production is still a conjecture. 

The only exclusive quarkonium production process that has been observed at B factories 
is e + e~ — > J ftp + t] c [3]. Right after the first observation, it was revealed that the measured 
production rate [3] is greater than the leading-order (LO) NRQCD results [4, 5] by an order 
of magnitude. Although the subsequent measured values have decreased [6, 7] in comparison 
with the first measurement, the discrepancy still remained. Through extensive theoretical 
studies on the corrections at the next-to-leading order (NLO) in the strong coupling a s 
[8] and the resummed relativistic corrections [9, 10], it was found that the interplay of 
the large QCD and relativistic corrections fill the gap between the theory and data [11- 
13] within uncertainties. A recent result for the order-a s t> 2 corrections also supports this 
conclusion [14]. However, some subtle issues that should be clarified still exist: The current 
experimental data for e + e~ — > J/ip + i] c contain the events with a fi + fi~ pair plus at least 
two charged tracks as the decay product of the r] c . If we take into account the t] c decay 
modes without charged tracks, then the actual cross section can be larger than the current 
data. While the size of the relativistic corrections has been estimated rather precisely with 
the aid of the resummation technique, the QCD corrections were computed only at order 
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a s . The amount of higher-order (HO) QCD corrections are not estimated yet. 

Another interesting issue is regarding the color-singlet LDMEs. The color-singlet LDMEs 
are usually determined by comparing experimental data with theoretical predictions for the 
electromagnetic decays of the S'-wave charmonia, J/ip — > / + /~ and t] c — > 77. It is known that 
the relativistic corrections to the charmonium decay rates are significant. In order to deter- 
mine the color-singlet LDMEs including the relativistic corrections to these electromagnetic 
decays, we need at least two constraints for each decay rate. One is the experimental decay 
rate. As for the other constraint, one may choose the generalized Gremm-Kapustin rela- 
tion [10] that relates the LDMEs of LO and HO [15-18]. As an alternative, one can use the 
measured rate T[J/ip — > LH] for the J/ip decay into light hadrons (LH) [13]. However, the 
numerical values for the LDMEs have a strong dependence on the choice: The LO LDME 
for J/ip at m c = 1.4 GeV is 0.440 GeV 3 [15] and 0.573 GeV 3 [13]. The corresponding values 
for v 2 are 0.225 [15] and 0.089 [13]. Therefore, the phenomenological study on the exclusive 
quarkonium production rate in e + e~ annihilation may provide us with a good chance of 
improving the accuracies in the determination of the color-singlet LDMEs. An accurate de- 
termination of the color-singlet LDMEs may lead to improving the determinations of various 
color-octet LDMEs that are involved in the inclusive productions of J/ip [19]. 

One may also consider exclusive S'-wave quarkonium processes in e + e~ annihilation other 
than the observed one e + e~ —> J/ip+i] c . In Refs. [20, 21], the process e + e _ — > J/tp+J/ip was 
first introduced as a possible source of contaminating the data samples for e + e~ — > J/tp + r] c . 
The angular-distribution analysis of the e + e~ — > J/ip + 7] c events at Belle showed that there 
were no such contaminations. In addition, the Belle Collaboration has not observed the 
signals for e + e~ — > J/ip + J/ip yet and only an upper bound of the cross section has been 
reported: a[e + e~ J ftp + J/ij)\ x B >2 [J/^} < 9.1 fb [6]. According to the QCD NLO 
correction in Ref. [22], the NLO corrected cross section is significantly smaller than the 
LO prediction. As is stated in Refs. [20, 21], the process is dominated by the photon 
fragmentation. The relativistic and QCD corrections to this contribution can be summed 
effectively by making use of the vector-meson-dominance (VMD) approach [23]. The cross 
section within the VMD approach is a[e + e~ — > J/tp + J/tp] = 1.69 ± 0.35 fb at the center- 
of-momentum (CM) energy a/s = 10.58 GeV [23], which is significantly smaller than the LO 
prediction and is consistent with the nonobservation of the process at Belle. The production 
of a charge-conjugation parity +1 quarkonium associated with a photon, e + e _ — > H+j, was 
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first suggested in Ref. [24] as a nice probe to the color-singlet mechanism of NRQCD and the 
convergence of relativist ic corrections, especially for rj c (2S). Later, both the QCD NLO and 
relativistic corrections to this process were calculated and found to be considerable [25, 26]. 

In this work, we investigate the resummed relativistic correction and its interplay with 
the QCD NLO correction to the exclusive processes e + e~ — > r] c + 7 and e + e~ — > J/ip + J/ip 
at y/s = 10.58 GeV. Because the QCD NLO and relativistic corrections to e + e~ — > J/ip + rj c 
are large, one may guess that the two processes listed above may also acquire significant 
corrections relative to the LO contribution. Our calculation reveals that the corrections 
are indeed large negative. The resummation of relativistic corrections enhances the cross 
section for e + e~ — >■ ^ c + 7 in comparison with the order-?; 2 relativistic correction. In the case 
of e + e~ — > J/ip + J/4>, the QCD NLO and relativistic corrections resummed to all orders 
in v make the cross section negative within the standard NRQCD factorization approach. 
By employing the VMD approach to compute the photon-fragmentation contribution, we 
cure the problem. The result shows that the QCD HO corrections may have significant 
contributions to e + e~ — > J/ip + J/ip. Therefore, these two exclusive processes may provide 
us with independent phenomenological constraints to the color-singlet NRQCD LDMEs for 
the S'-wave charmonia. 

The remainder of this paper is organized as follows. In Sec. II, we briefly describe the 
strategy of resumming relativistic corrections in a quarkonium process. In Sec. Ill, we com- 
pute the cross section for e + e~ — > t] c + 7, in which the QCD NLO and resummed relativistic 
corrections, and their interference are included within the standard NRQCD factorization 
approach. The corresponding cross section for e + e~ — > J/ip + J/ijj is given in Sec. IV. In 
the case of e + e~ — >■ J/ip + J/ip, that is dominated by the photon fragmentation, the result 
within the standard NRQCD approach gives negative cross section. As an alternative, we 
also employ the VMD approach to find an improved prediction. Finally, we conclude in 
Sec. V. 

II. STRATEGY OF RESUMMING RELATIVISTIC CORRECTIONS 

In this section, we list the NRQCD factorization formula for the cross sections of e + e~ — > 
r] c + 7 and e + e~ — > J/ip + J/ip, and summarize the strategy of resumming relativistic 
corrections. 
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A. NRQCD factorization formula 

According to the NRQCD factorization formalism [1], the cross sections for the exclusive 
processes e + e~ — > r\ c + 7 and e + e~ — > J/ip + J/ip can be expressed as: 



where (O^ ) is the NRQCD LDME for production of a quarkonium H from a cc n pair with a 
spectroscopic state n, F n and F m ^ n are dimensionless perturbative short- distance coefficients 
that are independent of the long-distance nature of H, and m c is the charm-quark mass. 

In general, the quantum number for cc n does not have to be the same as that for the 
quarkonium H and the numerical importance of a cc n channel relative to the color-singlet 
channel is classified in powers of v in the velocity-scaling rules of NRQCD. In the exclusive 
process (1), cc n channels with the quantum number identical to that of H contribute domi- 
nantly. Thus the series (1) can be well approximated by the relativistic corrections and the 
index n can be understood as the power in v 2 relative to the LO color-singlet NRQCD LDME 
(Oq). The expression (1) can further be simplified once we apply the vacuum-saturation 
approximation to express the production LDME in terms of the corresponding decay LDME 
{Ou)h'- {@n) ~ (2^+ l)(£*n)#, where J is the total-angular-momentum quantum number 
of H. In the Coulomb gauge, the relative-order- 1> 2n decay LDMEs for H = J/ip and i] c are 
expressed as 



(O n ) m = (JMA)|^(-f£) Ve*(A)x|0}<0| X t(_±£) <r.e(A)^|JMA)>, 

(2a) 



where n = a + b, 1 ^ and x are two-component Pauli spinor fields that create a heavy quark 
and a heavy antiquark, respectively, a % is a Pauli matrix, D is the spatial component of the 
covariant derivative, and A and e(A) are the helicity and polarization vector of J/ip, respec- 
tively. 2 The state \H) is normalized nonrelativistically: (H(P)\H(P')} = (2tt) 3 5 (3) (P-P'). 

1 If a b, the right side of Eq. (2) is understood to be the average with its Hermitian conjugate. 

2 In general, there arc additional contributions that depend on chromoelectromagnetic field operators that 




(lb) 



(la) 




(2b) 
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The LO color-singlet LDMEs are given by 

(O ) W = |(0| X V-6(A)^|J/^(A))| 2 , (3a) 
(Oo) Vc = \(0\x^\Vc)\ 2 ■ (3b) 

B. Resummation of relativistic corrections 

The amplitude of an S'-wave charmonium H produced associated with a particle a in an 
e + e~ annihilation can be expressed as 



n=0 L V ^ / 



M(q 2 



(q 2n ) H (O )f, (4) 



where M.(q 2 ) is the corresponding parton-level amplitude with the standard normalization 
in the NRQCD factorization approach in which the angular dependence of half the relative 
momentum q of the c and c has been averaged in the cc rest frame. {q 2n )H is the ratio of 
the HO LDME of order v 2n relative to the LO one. For H = J/i/j, 

, 2n , (0\xH-^) 2n <7-em\J/m) „ 

{q )m (0|xtcr.e(A)^|J^(A)> " W 

Eventually in the factorization formula (1), the long-distance factors in Eq. (4) are expressed 
in terms of the LDMEs in Eq. (2) that are, in principle, independent of the LO LDMEs (Oq)h 
in Eq. (3). In Ref. [10], the Cornell potential model was used to express the ratio {q 2n )n in 
terms of (q 2 )ii'- 

(q 2n )H = (q 2 ) n H , (6) 

that we call the generalized Gremm-Kapustin relation 3 . The relation (6) neglects the gauge- 
field contribution to the covariant derivative in the Coulomb gauge and neglects the spin- 
flipping interactions so that {q 2n ) Vc ~ (q 2n )j/^- Various applications and detailed descrip- 
tions of the uncertainties of applying the relation can be found in Refs. [10, 11, 15, 16, 27]. 



appear from relative order v A . At the relative order i> 4 such contribution can be expressed in terms of 
the operators of the same relative order in Eq. (2) [27, 28]. In this work, we neglect the gauge-field 
contribution in the Coulomb gauge and, therefore, the contributions of the chromoelectromagnctic field 
operators are also neglected. 
3 The original form of the Gremm-Kapustin relation is (q 2 )u ~ {n^H — 2m c )m c [29], where m# is the 
S'-wave quarkonium mass. 
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The generalized Gremm-Kapustin relation (6) allows one to resum a class of relativistic 
corrections to all orders in v. The resultant amplitude is 

A[e + e~ ^H + a}= M((q 2 ) H )(O )f . (7) 

It is straightforward to apply this resummation method to e + e~ — > r\ c + 7 and e + e~ — > 

j/i/> + j/ip. 

C. The VMD approach 

In Ref. [4], the importance of the photon- fragmentation contribution in the J/i/j pro- 
duction in e + e~ annihilation was considered and the competition between the QED and 
QCD contributions was first observed. The dominance of the photon fragmentation in 
e + e~ —> J/ij} + J/ip at the B factories was predicted in Refs. [20, 21, 23]. In spite of a sup- 
pression factor a 2 /oi 2 s relative to the QCD subprocesses of other exclusive two-quarkonium 
production in e + e~ annihilation, the photon-fragmentation contribution is largely enhanced 
by a kinematic factor [ v / s/(2m c )] 4 that is associated with the fragmentation of a photon 
into a cc pair. The rate is further enhanced by a factor of log[y/s / (2m c )] originated from 
the collinear emission of a virtual photon in the forward region. Similar enhancement of 
the photon-fragmentation contribution was predicted in inclusive J/ip production with suf- 
ficiently large transverse momentum px at hadron colliders [30, 31]. 

As is stated in Refs. [11, 23], the prediction of the photon- fragmentation contribution 
in e + e~ — > J/ip + J/ip can be improved by employing the VMD approach. We shall find 
later that after combining the QCD NLO and resummed relativistic corrections to e + e~ — > 
J/i/j + J/ip, the cross section becomes negative within the standard NRQCD factorization 
approach. As an alternative, we employ the VMD approach to cure this symptom. 

If the photon-fragmentation dominance is valid, then one can neglect the non- 
fragmentation contributions that make a gauge-invariant subset. In that case, the amplitude 
is factorized into the product of Ai[e + e~ — > 7* +7*] and the 7* — > J/tp fragmentation factor 
that involves the local electromagnetic current 

jr(x) = e c ec(x)Yc(x) (8) 

for a free cc pair created at a fixed point x. Here, e is the electron charge and e c is the 
fractional electric charge of the charm quark. In VMD, the long-distance factor for the 
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transition of the pair to J /if) is expressed as 



(J/^(A)|^(x = O)|O)=^/. 07 e^(A), 



(9) 



where gj/^-y is the effective J/^-photon coupling. Applying this vertex to J /if} — > e + e , one 
can determine the coupling gj/^ [11, 23] as 



If we can neglect the virtual-gluon corrections to the photon-fragmentation process that 
connect the two cc lines, then the vertex (9) collects the QCD HO and relativistic corrections 
effectively. We can make an educated guess that a virtual-gluon correction, that appears 
from relative order a 2 , is negligible because the contribution is suppressed by a 2 relative 
to the LO and the kinematic enhancement is reduced due to the stronger virtualites of the 
internal lines. Because we use the measured values for the decay rate Y\J/ip — > e + e~] and the 
J /if} mass mju, in evaluating the effective coupling (10), large theoretical uncertainties can 
be replaced with tiny experimental errors. We take the scale mj/^, which is the momentum 
scale at the vertex, in computing the fine structure constant a in Eq. (10). 

D. Interference 

In a previous analysis in Ref. [11] on e + e~ — > J/if> + r) c , the authors found that both QCD 
NLO and resummed relativistic corrections to the cross section were significant. It was also 
reported that the interference between the two corrections was about 26 % of the LO cross 
section. Thus we expect that the interference term may also be important in other processes 
that have large QCD and relativistic corrections. In addition to this interference there is 
another order-a s v 2 correction that can be computed directly from the order-a s diagrams 
by keeping the v dependence. Such order-a s t> 2 corrections were computed for J /if) — > 
[16], B c lv [32], r) c 77 and LH [18, 33], e+e" J/ip + rj c [14], h c , h b and rj b -> 
LH [34], and J /if} — > 37 [35]. Many of the corrections are tiny, for example, at most 0.3% 
in J /if) — > £ + £~ [16]. Some exceptions are h c — > LH [34] and J /if} — > 37 [35], which have 
significant order-a s t> 2 corrections. 

In this work, we improve the theoretical prediction for the production rates of e + e~ — > 
r/ c + 7 and e + e _ — > J /if} + J /if) by including the order-a s and order-?; 00 corrections, and 




(10) 
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FIG. 1: Feynman diagrams for e + e — > r] c + 7 at LO in a s . 

their interference. 4 The computation of the pure order-a s w 2 contribution, that is out of the 

scope of this work, is not included. The missing contribution is of the same order as the 

interference term that may bring in theoretical errors. Equation (42) of Ref. [11] includes 
the three contributions listed above to the total cross section <j tot : 

, 0"NLO — °"0 , ,js i\ m\ 

Ctot = &v°° + V "^ 00 F= = CT « 00 + V a v°°&0 [J^a 3 — J-J , [i-i-l 



where o"o, <7„oo, and ct^lo are the cross sections at the LO, order v°°, and order a s , re- 
spectively, and K aB = ctnlo/co- The approximation for the interference term in Eq. (11) 
is constructed under the assumption that K as and v°° corrections are independent of the 
polarization A of H . Violation of this assumption may cause additional theoretical errors. 



III. e+e — > T) c + 7 

In this section, we compute the cross section for e + e~ — > t] c + 7 at -y/i = 10.58 GeV by 
employing the strategies given in Sec. II. In e + e~ annihilation, this process proceeds with the 
parton process e + e~ —> cc( 1 S , q 1 ') + 7, where cc( 1 5'q 1 ') is the color-singlet spin-singlet S-wave 
cc pair. Because both QED and QCD have the charge- conjugation symmetry and r] c is an 
even charge-conjugation eigenstate, the process factors into the leptonic part e + e~ — > 7* and 
the hadronic part 7* — > cc( 1 S'q 1 ') + 7. At the LO in a SJ two Feynman diagrams contribute 
to this process as shown in Fig. 1. 

4 In the remainder of this paper, we use the notation v°° in order to indicate that a quantity is resummed 
to all orders in v. 
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A. Amplitude 



In order to construct the amplitude for e + e — > cc^S^) + 7, we first write the amplitude 
M. for the process e + (k 2 )e~(ki) — > c{p)c{p) + 7(P 7 ) as 

M = Wly* (12) 

where the momentum of each particle is written in the parentheses following each particle 
identity. Here, the leptonic current is defined by 

L» = -v(k 2 )Yu(h). (13) 

s 

The hadronic current is given by = m(c)C /j w(c) with = C^ u e* v ', e is the polarization 
vector of the photon, and the Dirac and color indices of C^ v are suppressed. The momenta 
p and p for the on-shell c and c, respectively, are linear combinations of the quarkonium 
momentum P and half the relative momentum q: 

p = ±P + q, (14a) 
p = \P-q. (14b) 

At the cc rest frame, P = [2E(q),0], q = (0,q), p = [E(q),q], and p = [E(q),—q], where 



E(q) = ^Jm 2 c + q 2 . 

The color-singlet spin-singlet contribution can be projected out from "H M by taking the 
trace after multiplying the projectors for the spin-singlet [IIi(p, j?)] and the color-singlet (tti) 
state. The corresponding projectors are defined by 

Hifap) = - 1 -(p-m c ) 1 5 if+2E(q)](p' + m c ), (15a) 

Ay/2E(q)[E(q) + m c \ 

m = -;Li, (15b) 

where N c = 3 is the number of colors and 1 is the unit color matrix. The spin projector 
(15a) contains the q dependence to all orders in v [28]. Then the cc^S^) contribution 
%^[cc{} S^)\ to the hadronic current is finally obtained by averaging over the directions of 
q in the cc rest frame: 

H,[cc{ l S^)\ = J ^ Tr [H^fad ® 7Ti] , (16) 
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where £l q is the solid angle of q and the trace is over the color and spinor indices. Now 
the complete q dependence in Eq. (16) is a function of q 2 . Carrying out the standard 
perturbative matching of the full-QCD amplitude for the hadronic current to the NRQCD 
counterpart and employing the resummation strategy stated in Sec. II B, we obtain the 
amplitude for e + e~ — > r\ c + 7 as 

A[e + e~ -> rj c + 7] = L 
B. Kinematics 



2N c ^E{q) 



n,[ccCs^)} 



(Oo)l /2 



(17) 



Because of the average over the angle for q that is required to project out the S*-wave 
state, additional consideration of the q dependence in the cc amplitude is needed. As is 
stated in Ref. [11], it is convenient to choose the e + e~ CM frame in which the final momenta 
are along the z-axis. In this CM frame, the momenta for e~, e + , r] c , and 7 are given by 

fs 

hi = -y (1, sin 0,0, cos 0) , (18a) 

k 2 = ^(l,-frin0,O,-cos0), ( 18b ) 
P = (E P , 0, 0, P C m) , (18c) 
P 7 = (Pcm, 0, 0, -P C m) , (18d) 

where 6 is the scattering angle, E P = [s + AE 2 (q)]/(2^), and P C m = [s — AE 2 (q)}/(2^). 
The explicit form of q in this CM frame is given by 

q= \q\ (jqPq cos9 q , sin 6^ cos <fi q , sin 9 q sin </> q , ^ q cos 9 q ) , (19) 

where |q| = \/ —q 2 and the two factors 7^ = Ep/[2E(q)] and (3 q = Pcu/Ep involve the 
boost from the cc rest frame to this e + e~ CM frame. 9 q and <j) q are the polar and azimuthal 
angles, respectively, of q* that is the explicit form of q in the cc rest frame: 

q* — \q\ (0, sin 9 q cos <p q , sin 9 q sin <p q , cos 9 q ) . (20) 



C. Cross section 

Summing over the final spins and averaging over the initial spins of the absolute square of 
the amplitude ^4[e + e~ — > ?7 c + 7] in Eq. (17), we obtain the spin-averaged squared amplitude 
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TABLE I: The cross sections gq, <t v °°, and <r t ot for e + e~ — > r] c ± 7 in units of fb at the scales 
fj, = m c and 2m c with the input parameters m c , {Oq)^, (q 2 ) Vc -> lli units of GeV, GeV 3 , and GeV 2 , 
respectively. The errors of the cross sections reflect the uncertainties of the NRQCD LDMEs only. 



m c 


(Oo) Vc 


<<z% 


00 


(TyOO 


CTtot 

[i = m c /i = 2m c 


1.4 


U-^-J' -0.105 


0.442 ±0.143 


00 + 21.1 

oo.z_ 20 


68.9™ 


c+13.6 c-r q+14.5 
0i -°-12.8 00 - -13.7 



|*4[e+e~ — 7- r] c ± 7]| 2 . Dividing it by the flux factor 2s, and integrating over the phase space, 
we find the cross section. 

In order to evaluate the cross section for e + e~ — > t] c ± 7 we first choose the numerical 
values for the input parameters. We choose a = 1/131 that is the running coupling constant 
evaluated at scale /1 = y/s [15, 24]. As is shown in the first three columns in Table I, we 
choose m c , (Oq) Vc , and (q 2 )n c from Refs. [15, 24]. The LO cross section o"o is consistent with 
previous results in Refs. [24, 25] 5 . We recall that a v oo is the order-a° contribution in which 
order-f°° corrections are resummed. The order-t>°° correction a v oc is about —17% of o"o. We 
have estimated the order-w 2 corrections from our resummed formula by taking the partial 
derivative with respect to (q 2 )n c numerically. Our result is — 75.6 {v 2 ) Vc , which is in good 
agreement with the value —75.5 (v 2 ) Vc obtained from the fixed-order calculation in Ref. [25]. 
We define K v % = a^ja^ and K v oc = <j vao /<j , where a v 2 is the cross section that includes the 
corrections of order a° s v 2 . The results are K v i = 0.80 and K v oc = 0.83 which imply that the 
contributions of order v A or higher enhance the cross section by about 4 %. This is a good 
signal that the v expansion converges very well. 

The factor K aa defined by ctnlo/^o is 0.77 for /i = m c and 0.82 for ji = 2m c , respectively 
[25]. Taking into account both the NLO corrections in a s and the relativistic corrections, and 
the interference between them, we find the total cross section a tot = 51.51^ 8 fb for )U = m c 
and <7tot = 55.3^^3 7 fb for yU = 2m c , respectively. For a luminosity of 1 (10) ab _1 at B (super 
B) factories, the expected number of events is about 5 x 10 4 (5 x 10 5 ). Since the branching 
fraction of the rj c decay into two photons is (1.78 ±0.16) x 10~ 4 [36], e + e~ — > ?7 C ±7 followed 
by rjc — > 77 might be observed at Belle II or super B factory, but the events might not be 

5 The errors of the cross sections reflect the uncertainties of the NRQCD LDMEs only. 
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triggered because no charged particles exist in the detector. Instead, the rj c — > KKir mode, 
whose branching fraction is (7.2 ± 0.6) % [36], may be useful to observe the e + e~ — > r\ c + 7 
events. The analysis of the photon energy spectrum in e + e~ — > KKir + 7 will be useful to 
detect this channel [24]. 

IV. e+e" J/V> + J/tp 

In this section, we consider the e + e~ — > J/ip + J/ip process at y/s = 10.58 GeV. This 
process proceeds through e + e~ annihilation into two virtual photons because the charge- 
conjugation parity of J/ip is — 1. At the LO in a s , there are four Feynman diagrams as shown 
in Fig. 2. The first two diagrams are called the photon-fragmentation diagrams, where each 
virtual photon evolves into a J/ip. The other two diagrams in the second row are called the 
nonfragmentation diagrams, where two virtual photons evolve into two J/tp mesons. The 
invariant mass of the virtual photon in the photon-fragmentation diagrams is fixed as mj/^, 
but the typical momentum scale of the virtual photons in the nonfragmentation diagrams 
is of order a/s/2. 

A. Amplitude of e+e" ->• cc( 3 s[ 1] ) + cc( 3 s[ 1] ) 

In order to calculate the relativistic corrections to the double J/ip production in e + e~ 
collisions, we first construct the amplitude for the production of two color-singlet spin-triplet 
S'-wave cc pairs in e + e collisions. The amplitude for e + {k,2)e (fci) — > c(pi)c(pi) + 
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c(p 2 )c(p 2 ) can be written as 

4 

M = Y,L^A^[ciCi + c 2 c 2 l (21) 
i=i 

where and are the leptonic and hadronic tensors for each diagram. Eventually, 
each Ci(pi)ci(pi) pair evolves into a J/ift meson. 

The amplitude A^ U [P\, q\, Xi, P 2 , q 2 , A 2 ] for the production of two cc( 3 Sf^) pairs from the 
two virtual photons can be obtained by projecting the amplitude A^ v \c\Cx + c 2 c 2 ] into the 
color-singlet spin-triplet channel, where Pi, qi, and Aj are the four- momentum of the CjC; 
pair, half the relative momentum of q and q, and the polarization vector of the QCj pair, 
respectively, with the following relations: 

Pi 

Pi = Y + qu ( 22a ) 

Pi = y - Qi- (22b) 
The spin-triplet projector for qq pair, which is valid for v°° corrections, is given by [28] 

mPi.PiAi) = - nr „ / 77^ : , r($- ^ [/> H- 2E{*)] ™e) , (23) 

4V2E(^)[E(^) + m c J 

where ej(Aj) is the polarization vector of the qq pair with the polarization Aj and -E^i) is 
the energy of the charm quark or antiquark in the QCj rest frame. We note that the CjCj pair 
in A^ U [P\, qi, Ai; P 2 , q 2 , A 2 ] is not a pure S'-wave state. 

The S'-wave contribution can be extracted by averaging the amplitude over the polar and 
azimuthal angles of qi, which is the spatial vector of qi, in the CjCj rest frame: 

^[(cici)( 3 sl 1] ) + (c 2 c 2 )( 3 sf)] ee J ^^A^[P 1 ,q 1 ,\ 1 ;P 2 ,q 2 ,\ 2 }. (24) 

Now the qi dependence of the quantity in Eq. (24) is a function of q\ and q\. Here, Qi is 
the solid angle of qi in the CjCj rest frame. As we did in the previous section, we do not find 
the short-distance coefficients by expanding Eq. (24) in terms of qf. Instead, we employ 
the formula (7) to obtain the amplitude in which relativistic corrections are resummed: 



AUJ/1> + J/1>] 



— 7 =^_=^[(c 1 c 1 )( 3 M 1] ) + (c 2 c 2 )( 3 sf)] 



<?I =92 = (<J 2 ).//</, 

x(Oo)^(Oo)^- (25) 
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B. Kinematics 

We first compute the cross section for e + e~ — >■ J/ip + J/ip within the standard NRQCD 
factorization approach. Then, as an alternative, we employ the VMD approach for the 
photon-fragmentation diagrams while keeping the standard NRQCD factorization approach 
for the nonfragmentation diagrams as in Ref. [23]. Note that we use the expression Pi 
for the CjCj four-momentum depending on the approach: In the NRQCD approach, we set 
Pf = [2E(qi)] 2 , while we set Pf = rn 2 ,, in the VMD approach. 



1. Kinematics for the NRQCD factorization approach 

In the calculation within the standard NRQCD factorization approach, the four-momenta 
of e and CjCj pairs at the e + e~ CM frame are given by 

(26a) 

(26b) 
(26c) 
(26d) 

<(<h) 2 + 4E(q 2 ) 2 ]/(2^). and 
\P\ = A 1/2 (s, AE( qi ) 2 , 4E(q 2 ) 2 )/(2^) with A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2bc - 2ca. The 
polarization vectors of the CjCj pairs are given by 

el(0) = * t (\PlOAEi), (27a) 



h = 


^(l,sin0,O,cos0), 


h = 


—— (1, — sin^, 0, — cosi 


Pi = 


(^i,0,0,|P|), 


P 2 = 


(E 2 ,0,0,-|P|), 


AE(q 2 


) 2 ]/(2^), E 2 = [s- 



E 2 - \P\ 2 
1 

VE'i - \p\ 2 



e*(0) = —=L==(-\P\,0,0,E 2 ), (27b) 



el(±) = - e *(=F)==F^(0,l,=Fi,0). (27c) 

In each rest frame of the qq pair, half the relative momentum q* is given by 

Qi = \Qi\ {0,sm9* coscj)*,sm9* smcf)*,cos9*) , (28) 

where 9* and <fi* are the polar and azimuthal angles of q*, respectively, and \qi\ is the 
magnitude of q* in the c,Cj rest frame. Boosting Eq. (28) from the rest frame to the 
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e + e CM frame, one obtains 



Qi — \Qi\ (liPi cos s i n $1 cos0^, sin 6{ sin <f>l, 71 cos 6{) , (29a) 
<72 = | (72 1 (T2 cos ^2 ) sin 02 cos 02 ; sin 9 2 sin 02 ; T2 cos ^2 ) 5 (29b) 

where 7i = Ei/[2E{ qi )], ft = |P|/£i, and (5 2 = -\P\/E 2 . 



2. Kinematics for the VMD approach 

In the calculation within the VMD approach, the four-momenta and polarization vectors 
for the CjCj pairs at the e + e~ CM frame are given by 

P* = (£f,0,0,|P fr |), (30a) 
P 2 fr = 0,0,-|P fr |), (30b) 
e f» (0 ) = 1 (|P fr |,0,0,Ef), (30c) 



_ |p 

e f»(0) = 1 (-1^1,0, 0,JSf), (30d) 

yjEf - |P fr | 2 

ef(±) = -ef (=p) = T -L(0, 1,=F*,0), (30e) 

where Ef = E% = and |P fr | = (s — 4m 2 7 ^) 1 / 2 /2. Here, we use the superscript fr to 

denote that the variables are for the VMD approach. The four-momenta for the e ± are the 
same as those in Eqs. (26a) and (26b): kf = fcj. Because the photon- J/ip vertex (10) already 
contains the relativistic corrections to the photon-fragmentation diagrams, we do not take 
the average over the angles of q± and q 2 . 
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C. Numerical results for NRQCD prediction 

Applying the method described in Sec. IV Bl, we carry out the numerical computation 
of the cross section a[e + e~ —> J/ip + J/ip] at the B factories within the standard NRQCD 
factorization approach. We take a = 1/131 and vary the charm-quark mass as m c = 1.4 
and 1.5 GeV. The results are summarized in Table II depending on the input parameters. 
The LO NRQCD LDME (Oo)jm and the ratio (q 2 )j/^ are listed on the second and third 
columns of Table II, respectively. The values for (O )jm and (q 2 )j/^p at m c = 1.4 GeV 
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TABLE II: The cross sections ctq, cr v °o, and <7tot f° r e + e~ —> J/ip + J /if) within the standard 
NRQCD factorization approach in units of fb at the scales [i = m c , 2m c , and y/s/2 with the input 
parameters m c , (Oo)jm, {q 2 )j/,p, in units of GeV, GeV 3 , and GeV 2 , respectively. The errors of 
the cross sections reflect the uncertainties of the NRQCD LDMEs only. 





{°o)j/il> 




o"0 


(J v oo 


H = m c 


O"tot 

fj, = 2m c 


H = yfi/2 


1.4 


n 44n+ - 067 

U.44U_ 055 


0.441 ±0.140 


13.39±1;P 


5-56±l;II 


-12.74±l;iS 


9 co+0.56 


-0.89^ 


1.5 


U.^CiD_ 054 


0.442 ±0.140 


8.09iJ:2 


3.66l o: | 




-1.36 ±0.33 


-0.43l°;i 



are taken from Ref. [15]. The values at m c — 1.5 GeV are obtained by applying the same 
method described in Ref. [15]. 

The result for the LO cross section o"o is listed on the third column of Table II. Our 
NRQCD prediction for cto is greater than a previous result o"o = 6.65 ± 3.02 fb in Ref. [20]. 
This is mainly because our value for (Oo)jm is greater than that (0.335 ± 0.024 GeV 3 ) of 
Ref. [20]. Our short-distance coefficient at LO in a s agrees with those in Refs. [20, 22] and 
the numerical results also agree once we use the same numerical value for the LDME. 

The cross section a v <x> in which the relativistic corrections are resummed to all orders 
in v is listed on the fourth column of Table II. The resummed result (JyOO IS SI gnificantly 
smaller than the LO prediction <jq leading to K v oo = 0.42 (0.45) at m c = 1.4 (1.5) GeV. As 
we have mentioned, we can extract the order-?; 2 relativistic corrections from the resummed 
formula by varying {q 2 )j/^, for each J/ip numerically. The correction is — 28.4 (v 2 ) j/^ fb 
( — 17.4 (v 2 )j/^fb) at m c = 1.4 (1.5) GeV for each J/tf>, where (v 2 )j/^, = (q 2 ) j/^/m 2 c . Taking 
into account the fact that there are two J/ip mesons, we find that K v 2 = 0.05 (0.15) at m c = 
1.4 (1.5) GeV. While the order-f 2 corrections are large negative, the relativistic corrections 
of order v 4 or higher significantly enhance the cross section. 

The authors of Ref. [22] computed the cross section <tnlo 

for e+e" -> J/ip + J/ip at NLO 
in a s , which has a strong dependence on m c and the renormalization scale \i. According to 
the results, the QCD NLO corrections are large negative, leading to negative cross sections 
at fi = m c for both m c = 1.4 and 1.5 GeV. We can extract the factor K as = onlo/^o from 
the results in Table I of Ref. [22]: At m c = 1.4 (1.5) GeV, K as = -0.367 (-0.314), 0.057 
(0.077), and 0.253 (0.248) for /x = m c , 2m c , and y/s/2, respectively. 
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Here, we notice that both the QCD NLO and order-?; 2 relativistic corrections are large 
negative and significantly reduce the cross section for e + e" — > J/'if) + J /'if). Thus, the 
interference between them must be large positive. In addition, the relativistic corrections of 
order v A or higher is large positive. Therefore, we expect that the interference contribution 
can be significant. 

A rigorous computation of the interference contribution requires the information of the 
factor K aa for each polarization of two J /if) mesons, which is not available yet. Therefore, we 
ignore the dependence on the J/if> polarization and use the formula (11) for <r tot to include 
the QCD NLO corrections, relativistic corrections, and their interference. Our results for 
0"tot depending on the input parameters are given on the last three columns in Table II. 6 
Every entry for <7 to t is negative, which is unphysical. We expect that one can resolve this 
problem only after including the corrections of orders a 2 , a s v 2 , or higher. 

D. Numerical results with the VMD approach 

In this section, we provide the numerical results for the cross section for e + e~ — > J/ip+J/ip 
in which the photon-fragmentation contribution is computed with the VMD approach. Ac- 
cording to the NRQCD predictions in Refs. [20, 21], the nonfragmentation contribution 
to e + e~ — > J /if) + J /if) cross section is 0.9%, and the interference between the photon- 
fragmentation and nonfragmentation contributions is —15.4 %. This reveals that the photon- 
fragmentation contribution dominates in the process. In that limit, the VMD approach 
provides a good approximation for the photon- J /if) vertex in the photon-fragmentation dia- 
grams. The coupling in Eq. (10) collects both the QCD NLO and relativistic corrections to 
the photon- J /if) vertex effectively and, therefore, the theoretical uncertainties can be greatly 
reduced [11]. Note that we cannot employ this method to compute the nonfragmentation 
diagrams and we apply the conventional NRQCD approach that is described in the previous 
section. This hybrid strategy can be justified because each contribution makes a separate 
gauge- invariant subset. 

In estimating the coupling gj/^ 7 in Eq. (10), we take T[J/ip — > e + e~] = (5.55 ± 0.14 ± 
0.02) keV [36] and a(fi = mj/^) = 1/132.6, where the scale [i corresponds to the momentum 

6 At /i = m c , K aa is negative and, therefore, Eq. (11) is inapplicable to computing o- t ot- In this case, we 
take o-tot = cr„oc + {K as - l)a . 
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TABLE III: The cross sections , cr^, and cr^ ot for e + e~ — >■ J/ip + J/ip in units of fb, in which the 
photon-fragmentation contribution is computed with the VMD approach when gj/^y = 0.832 ± 
0.010 GeV 2 . The description on the input parameters is the same as that in Table II. 





{Oo)j/ 1 p 


(q 2 )j/i, 




fT fr 


[i = m c 


^ = 2m c 


M = v^/2 


1.4 


U-^HJ-0.055 


0.441 ±0.140 


1.81 ±0.13 


1.80 ±0.13 


1 n5 +0 - 22 
i.uo_ 025 


Q9+0.18 

i.oz_ 20 


1 qq+0.17 


1.5 


U.^OD_ 054 


0.442 ±0.140 


1.85 ±0.13 


1.84 ±0.13 


1.35 ±0.20 


1 4Q+0.17 
-■■•^-0.18 


1.56 ±0.16 



transfer at the vertex. In other vertices, we take a = 1/131. The resultant numerical value 
is gj/ ^ = 0.832 ± 0.010 GeV 2 . 

Our results for the cross sections (Tq, cr^, and cr£, t are listed in Table III. Here, the 
superscript fr indicates that the photon-fragmentation contribution is computed with the 
VMD approach while the nonfragmentation diagrams are computed in the conventional 
NRQCD factorization approach. The description on the input parameters is the same as 
that in Table II. The uncertainties of the cross sections are dominated by the experimental 
errors in T[J/ip — > e + e~] through gj/^. The LO cross sections <Jq at m c = 1.4 and 1.5 GeV 
agree with a previous result 1.69 ± 0.35 fb in Ref. [23] within errors. Sources of this small 
difference are as follows: While we use 2Ei(qi) for the invariant mass of the cc pair for the 
nonfragmentation diagram, the authors of Ref. [23] used the physical mass mj^ to normalize 
the heavy-quarkonium state. Another source of the difference is that we use a = 1/132.6 to 
determine gj/^, whose scale is taken to be mj/^, but in Ref. [23], a = 1/137 was used. 

We notice that the LO cross section in Table III is significantly smaller than o"o shown 
in Table II. This is because the relativistic corrections to all orders in v are already included 
in the fragmentation contribution to that was computed with the VMD approach. In 
addition, and a^ 2 are essentially the same as <Jq. As a result, K% « K% w 1. This 
indicates that the relativistic correction to the contribution not coming from photon frag- 
mentation is negligible. 

It is important to notice that, in the VMD approach, the QCD and relativistic corrections 
to the photon- J/ip vertex have been effectively resummed to all orders. Therefore, one must 
be careful to avoid double counting of the QCD NLO corrections when one combines the 
NRQCD one-loop corrections and the relativistic corrections. In each photon- J/ip vertex of 
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the photon-fragmentation diagrams, there is a QCD NLO correction factor [1 — 16a s /(37r)] 
at the cross-section level. Therefore, in order to subtract the amount that is contained in 
the NRQCD one-loop corrections in Ref. [22], we add 32a s / (3ir) to the factor K ag extracted 
from the result in Ref. [22]. The resultant factor is very close to 1: K£ = 0.944 (0.939), 
0.964 (0.956), and 0.969 (0.964) for /i = m c , 2m c , and y/s/2, respectively, at m c = 1.4 
(1.5) GeV. 

Because most of the QCD and relativistic corrections are already contained in a^, we 
have found that ~ 1 and m 1. In this limit, the interference between the QCD 
and relativistic corrections is negligible. Therefore, we take a^ ot « K^ctq + (K£ — l)ao- 
The results for cr[ ot are listed on the last three columns of Table III. It is remarkable that 
every entry for a^ ot listed in Table III is positive in contrast to the corresponding values 
in Table II. This is a strong indication that the corrections of order higher than a s v° may 
resolve the problem that the NRQCD prediction gives negative cross section. In addition, 
of ot is insensitive to the variation of the scale fi from m c to y/s/2. Our prediction for 
cr[ ot is consistent with the nonobservation of the process at Belle with an upper bound 
a[e + e~ -> J/i/j + J/ip] x B >2 [J/^} < 9.1 fb [6]. 

V. CONCLUSION 

We have computed the cross sections for the two exclusive production processes, e + e~ — > 
f] c + l and e + e~ — > J/ip + J/ip at = 10.58 GeV, in which a class of relativistic corrections 
are resummed to all orders in v. By including the available QCD corrections of order a s 
and the interference between the QCD and resummed relativistic corrections, we further 
improved the corresponding theoretical predictions. 

The NRQCD prediction for the cross section of e + e~ — >■ rj c + 7 is listed in Table I. The 
value <7tot ~ 50 fb, that is about 60-70% of the LO cross section ctq, contains the QCD 
NLO and resummed relativistic corrections, and their interference. In this process, both the 
QCD NLO and order-?; 2 relativistic corrections are negative. The relativistic corrections of 
order-w 4 or higher are negligible. The value for cr tot is still large enough to be measured at 
B factories or super B factories, for example, through the t] c — > KKn decay mode. 

In the case of e + e~ — > J/ip + J/ift, which is dominated by the photon fragmentation, 
we have found that the standard NRQCD prediction for crtot, shown in Table II, is nega- 
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tive for every choice of input parameters. As an alternative, we have employed the VMD 
approach to include the QCD and relativistic corrections to the photon- J/ip vertex in the 
photon-fragmentation diagrams. The VMD result a^ ot ~ lfb, given in Table III, indeed 
resolves the problem of negative cross section and the prediction is insensitive to the scale [i. 
This is a strong indication that corrections of order higher than a s v may have significant 
contributions. The predicted cross section for e + e~ — > J /if) + J /if) is consistent with the 
nonobservation of the events at Belle. 
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